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This paper provides a review and fundamental physical interpretation for the effective densities 
and moduli of acoustic metamaterials. We introduce the terminology of hidden force and hidden 
source of volume: the effective density or modulus is negative when the hidden force or source of vol¬ 
ume is larger than, and operates in antiphase to, respectively, the force or volume change that would 
be obtained in their absence. We demonstrate this ansatz for some established acoustic metamateri¬ 
als with elements based on membranes, Helmholtz resonators, springs and masses. The hidden force 
for membrane-based acoustic metamaterials, for instance, is the force from the membrane tension. 
The hidden source for a Helmholtz-resonator-based metamaterial is the extra air volume injected 
from the resonator cavity. We also explain the analogous concepts for pure mass-and-spring systems, 
in which case hidden forces can arise from masses and springs fixed inside other masses, whereas 
hidden sources—more aptly termed hidden expanders of displacement in this case—can arise from 
light rigid trusses coupled to extra degrees of freedom for mechanical motion such as the case of 
coupling to masses that move at right angles to the wave-propagation direction. This overall picture 
provides a powerful tool for conceptual understanding and design of new acoustic metamaterials, 
and avoids common pitfalls involved in determining the effective parameters of such materials. 


Acoustic metamaterials are man-made structures de¬ 
signed to manipulate the propagation of sound in ways 
not available in naturally occurring materials. The un- 
derstandirm of negative constitutive parameters in such 
materials. |ll-l34| one of their most exotic features, has to 
date relied more on engineering-based concepts than on 
universal physical principles. This can lead to confusion 
in assigning effective parameters to a given system. Here 
we present a fundamental physical picture to account for 
the effective densities and moduli of acoustic metamateri¬ 
als that allows their intuitive yet precise understanding, 
and at the same time allows their unambiguous deter¬ 
mination. This picture is based on hidden forces and 
hidden sources of volume. The former, resulting in an 
effective mass or density, involve local, time-dependent 
non-apparent forces that provide a net force on the unit 
cell of the acoustic metamaterial. The latter, resulting in 
an effective elastic modulus, involve time-dependent hid¬ 
den sources of volume or displacement that only produce 
pairs of forces acting equally and oppositely on either 
side of the unit cell. 


We first illustrate our approach with examples of sys¬ 
tems exhibiting an effective mass based on membranes 
in tubes as well as on mass-and-spring models (the latter 
providing elements to model three-dimensional (3D) solid 
acoustic metamaterials). An example of a metamaterial 
based on membranes combined with masses and springs 
is also elucidated. We then provide examples of sys¬ 
tems exhibiting an effective modulus based on Helmholtz- 
resonators in tubes as well as based on mass-and-spring 
models combined with light rigid trusses coupled to ex¬ 
tra degrees of freedom for mechanical motion. Systems 
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exhibiting both effective density and modulus, includ¬ 
ing the possibility of double-negative parameter behav¬ 
ior, are also discussed. 


EFFECTIVE DENSITIES: HIDDEN FORCES 

In this section we introduce the concept of “hidden 
force” in order to understand effective mass in systems 
that contain non-apparent mechanical elements. [33 This 
approach is first explained by means of simple examples 
from mechanics. We then illustrate the concept of ef¬ 
fective density for several different acoustic metamateri¬ 
als involving unidirectional propagation: a membrane- 
based metamaterial, a solid-matrix metamaterial ap¬ 
proximated by spring-coupled masses, and also a new 
membrane-based metamaterial that also includes masses 
and springs. 


Simple mechanical systems 

Consider a wheel of mass M, radius R and moment 
of inertia / rolling in the x direction without slipping 
on a flat, horizontal surface, as shown in Fig. 1(a). A 
horizontal force F is applied to the axis of the wheel at 
its center. The application of Newton’s laws allows one 
to derive the acceleration: x = F/[M{\ -\-1/MR^]). The 
effective mass of wheel can be defined as 


Meff = F/x, (1) 

giving, in this case, M^ff = M{1 -\- I/MR^). Let us 
define the hidden force Fh to be the backwardly-directed 
frictional force on the wheel rim. The acceleration x is 
smaller than that expected for a non-rotating mass of 
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HIDDEN FORCE/EXTERNAL FORCE F^/f 


FIG. 1: (a) Showing external and hidden forces, F and Fh, 
for the case of a rolling wheel, (b) shows a normalized plot of 
the effective mass as a function of the hidden force, (c) shows 
F and Fh for a prototype of a system exhibiting negative 
effective mass: a simple harmonic oscillator. 


resonance frequency. Substituting into Eq. 

= (3) 

The displacement, obtained by integrating x = F/Meff, 
can be clearly seen to oscillate with large amplitude near 
resonance (at wq) because Mg// becomes very small. 
This ansatz describes all physical quantities correctly and 
quantitatively. Another example is the abrupt shift of the 
phase of the displacement by tt with respect to the driv¬ 
ing force as the frequency passes through the resonance. 
This can immediately be understood from Eq. m , since 
the sign of Mgff changes at w = wq. Negative a 

consequence oi Fh < —F, implies in the case of sinusoidal 
excitation that the magnitude of the hidden force is not 
only in antiphase with but also has a magnitude that is 
larger than that of the applied force. As the limit uj = 0 
is approached, the effective mass tends to —oo because 
the required force F for a given oscillation amplitude be¬ 
comes tiny in comparison with the oppositely directed 
spring force Fh- As the limit w = oo is approached, the 
effective mass tends to M because the (inertial) force F 
required for a given oscillation amplitude becomes very 
large in comparison with the spring force Fh- 


Membrane-based acoustic metamaterial 


the same magnitude, i.e. x < F/M. By introducing the 
effective mass as so defined, it is thus possible to obtain 
the correct acceleration from one simple equation. 

In general, if an external force F is applied to a mass 
M, then, owing to the specific mechanism involved, a 
hidden force Fh may also act on the mass (assumed here 
to be collinear with F). So the acceleration becomes 
X = {F + Fh)/M, and the effective mass is given by 


Meff 


M 

l + Fh/F' 


( 2 ) 


A plot of Mf.ff/M vs Fh/F is shown in Fig. 1(b). No¬ 
tably, as Fh approaches —E, M^^ff becomes infinitely 
large. Also, becomes negative when Fh < —F. 

A prototype of a system exhibiting negative M^f / is a 
simple harmonic oscillator consisting of a mass M at¬ 
tached to a rigid wall by a spring, as shown in Fig. 
1(c): the hidden-force picture starts by ignoring the 
presence of the spring and regarding the system as a 
free mass subject to a hidden force Fh = —kx, where 
k is the spring constant and x the displacement. In 
the case when the external driving force is sinusoidal 
at angular frequency lo, F = Eoexp(—iwt), the accel¬ 
eration can be calculated from the equation of motion, 
Mx = Fq ex-p{—iujt) + Fh- Using the harmonic expres¬ 
sion X = xq exp{—iujt), we then obtain the hidden force 
Fh = —kx = —FuJq/{uJq — w^), where luq = \/k/M is the 


Acoustic metamaterials, consisting of arrays of res¬ 
onators, naturally fit into this hidden-force picture. 
Take the example of a ID membrane-based acoustic 
metamaterial, S, S 11 > 0, 0, 0, 0| which supports 


wave propagation down its length, as shown schemati¬ 
cally in Fig. 2(a). It consists of a cylindrical air-filled 
tube containing taut membranes at regular intervals. 
Consider a particular unit cell. Its center of mass M is 
subject to two kinds of forces: the applied force F = S/S.p 
from the two adjacent cells, where S and Ap are the 
cross-sectional area of the tube and the pressure differ¬ 
ence across the unit cell. The hidden force is Fh = —kmF 
where km and ^ are the membrane spring constant |8| 
and the displacement of the unit-cell center of mass.|51j 
The equations governing oscillatory motion are the same 
as those governing the system of a mass connected by 
a spring to a rigid wall that we just treated: = 

SAp -|- Eft,, where M = pqSD -|- Mm is the mass of a unit 
cell. Here, po is the density of air, D is the unit-cell length 
and Mm is the mass of the membrane. For sinusoidal ex¬ 
citations, E = SAp = -I- km^ = — kmluj"^), so 

I = E/[M(1 - w^/w^)], where wq = s/km/M. 

According to Eq. ([1]), the unit-cell effective mass is 
Meff = M(1 — Wq/w^), which has exactly the same form 
as Mef f for a simple harmonic oscillator. This treatment, 
based on lumped-elements, obviously ignores vibrational 
resonances of the membranes higher than the fundamen¬ 
tal mode. 01 
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FIG. 2: (a) Schematic diagram of a ID membrane-based acoustic metamaterial, made up of periodically-spaced membranes in 
a tube containing air. Normalized plots (b), (c), (d) show the frequency as a function of the density, phase velocity and wave 
number, respectively. In (c) and (d) the solid and dashed lines refer to the cases of real and imaginary values for the phase 
velocity and wave number, respectively. 


For each unit cell indexed by j, the effective mass can 
be defined as 


Meff = - 


^3 


Sjpj -Pj+i) 

4 


(4) 


where Fj is the force acting on the left-hand side of the 
unit cell j, whereas is that acting on the right- 

hand side. (Alternatively, if we define fj as the force on 
the left-hand side and +fj+i as the force on the right, we 
obtain M^ff = {fj + which demonstrates more 

clearly that F = Fj — Fjj.i = fj + fj+i is the net force 
on unit cell j. However, in this paper we adopt the Fj 
notation because it is more convenient for the analysis 
of mass-and-spring models.) The effective density can be 
defined by pe// = where V = SD is the volume 

of the unit cell, which is a useful concept for the case 
Z? <C A, where A is the acoustic wavelength. It is this 
limit that applies to metamaterials, as opposed to the 
case for phononic crystals for which H ~ A. Therefore 


Peff = p' 



(5) 


where the average density p' is given by p' = M/V. 

In order to derive the system’s effective bulk modulus 
He//, consider the volume V between two adjacent mem¬ 
branes comprising parts of both the /th and (j — l)th unit 


cells. The non-equilibrium component of pressure in this 
volume, Vj .is related in general to the deviation AV) 
(which we define to refer to this same volume) from the 
equilibrium volume V as follows: 


AVi 

P3 — — -Beff- 


g fe - O-l) 

SD 


( 6 ) 


where Hg// is the effective modulus of the system. In 
the present geometry, only the (adiabatic) bulk modu¬ 
lus Bq of the air in the tube affects the pressure-volume 
relation. 


p, = (7) 

and so Hg// = Bq. Equation ([5]) represents continuity 
combined with the equation of state, and can be rewritten 
in the form 


lAV, 

Pj V Pj D 


( 8 ) 


In spite of the introduction of membranes in the tube 
with spring-like properties, they evidently do not change 
the effective modulus Bq of this acoustic metamaterial. 

One can also access the acoustic dispersion relation as 
follows. In the metamaterial limit, i.e. H A, the time 
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derivative of Eq. ([S]) leads, for an arbitrary point in the 
tube, to 


p = —B 


du 


(9) 


where u = ^ is the acoustic particle velocity and x is 
the distance along the tube. From the definition Pe// = 
Meff /(AD) and Eq. (U), we also have 

-( 10 ) 

The wave equation is easily obtained by combining Eqs. 
© and (fTUl) : 


■■ _ p d^P 
PeffP-Bog^^. 


( 11 ) 


Substituting p = po exp[i(gx — wt)], where q is the wave 
number, leads to the dispersion relation = uj^pef f /Bq. 
For the example in question. 


q = UJ 


Peff 


Bn 



1 - 


1/2 


( 12 ) 


A plot of the frequency dependence of the effective den¬ 
sity, and also plots for |up|, where Vp is the phase velocity 
(vp = uj/q) and jgl, where q is the wave number, are 
shown on normalized scales in Fig. 2(b)-(d). The phase 
velocity becomes infinite when the effective density van¬ 
ishes at a; = Wo- This situation is useful for applica- 

At 


tions in extraordinary acoustic transmission. [21 


lower frequencies, where Pe// is negative, the waves are 
damped [Im(q) > 0, dashed lines in (c) and (d)], whereas 
at higher frequencies, where the Pe// is positive, they are 
undamped [Im(g) = 0, solid lines in (c) and (d)]. The 
frequency region of damping, sometimes referred to as a 
metamaterial band gap, is useful in practice for applica¬ 
tions to the absorption of noise. 

One of the most curious features of these results is how 
membranes with spring-like properties only contribute to 
the effective density, and do not influence effective mod¬ 
ulus. We shall see later that this can also be explained 
by the fact that hidden sources of volume rather than 
hidden forces influence the effective modulus. 


Mass-and-spring analogy for a solid-matrix acoustic 
metamaterial 


shown in Fig. 3(b). The unit cell can be modelled by 
a core of mass Me surrounded by a hollow mass M, 
with an internal connection between the masses made 
up of collinear springs of constant fc/2. Springs of con¬ 
stant fco connect the masses M externally. The external 
springs transmit applied forces F to a particular unit 
cell, where F = Fq exp(—iujt), whereas the hidden force 
Fh = —k(^ — rf) acts on mass M through the inter¬ 
nal springs, where ^ and p are the displacements of the 
masses M and Me, respectively. Equations of motion for 
these parameters are M^ + k(^ — rj) = Fq exp(—iwt) and 
Mep + k(r] — ^) = 0 respectively. This clearly shows the 
origin of the hidden force as the vibration of the mass 
Me- Solving the coupled equations for sinusoidal motion, 
we obtain the effective mass of a unit cell in the form 


Meff = F/i =M + 


Me 


1 — w^/wj 


= M 1 + 


Wn — (JJ 


(13) 


where wi = wo\/l -I- Me/M and wq = yJk/Me (distinct 
from Wq in the above membrane problem). A plot of 
the normalized value of Meff (i.e. normalized effective 
density) vs frequency is shown in Fig. 3(c), for the case 
wi/wo = 2. One can see that as the frequency increases 
through w = Wo, there is a transition from inhnitely pos¬ 
itive to infinitely negative Meff. Zero Meff occurs at 


W = Wi. 

The effective modulus can be derived in a similar way 
to that for the membrane-based metamaterial. By anal¬ 
ogy with Eq. ©, 


= -keffi^j - O-i) = -D^Eeff ^^ (14) 

where keff is an effective spring constant and Eeff the 
effective Young’s modulus of the system. Alternatively, 

Eeff~~~^(ij~^j-^)- ( 1 ®) 

For the present ID mechanical model, the compressive 
force Fj is provided by the spring of constant fcg, i-®- Fj = 
-koi^j so keff = ko and Eeff = kn/D. Clearly, 

the effective modulus of this structure not affected by 
the internal structure of the mass. As in the case of the 
membrane-based metamaterial, the effective modulus is 
positive and frequency independent. 

The acoustic dispersion relation can be derived by ap¬ 
plication of Newton’s second law for sinusoidal variations: 


Now consider a metamaterial composed of a cubic ar¬ 
ray of mechanical resonators embedded in a compliant 
solid matrix, as illustrated by the unit cell in Fig. 3(a). 
Such a system consisting of rubber-coated lead balls was 
reported to exhibit strong sonic transmission loss due 
to negative density, [ij Consider a ID mass-and-spring 
analog of such a matrix in the form of a chain. |40l| as 


d'^F 

Mefftj = ko(^j-i - ^j) - koi^j - ^j+i) W koD^-^, 

""(16) 

where we have assumed, as before, that D X. Making 
use of Peff = Meff/D^ and Eeff = kn/D, we may write 

Pf^ff^ — Eeff ’ 


(17) 
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(a) 


(b) ^ 




UNIT CELL 



- D -► X 

UNIT CELL] 




PHASE VELOCITY |Vp|/(ko/p'D) WAVE NUMBER |q|/((BoP'D/ko) 


FIG. 3: (a) Unit cell of a solid-matrix acoustic metamaterial, modelled by a spherical core of mass Me surrounded by a separate, 
concentric rigid shell of mass M. (b) mass-and-spring analogy of a ID chain of unit cells. Fj is the compressional force in the 
spring of the left-hand side of unit cell j. Normalized plots (c), (d), (e) show the frequency as a function of the density, phase 
velocity and wave number, respectively. In (d) and (e) the solid and dashed lines refer to the cases of real and imaginary values, 
respectively, for the phase velocity and wave number. Normalized parameter uji/ujq = 2 is chosen for these plots. 


where 


Peff = P 


M 

M + Me 



(18) 


The constant p' = (M -|- Me)/D^ is the average density 
of the matrix. Substituting ^ exp[i(ga: — cut)] leads 
to the dispersion relation q = uj^peff/Eeff, or 


q = 



2 2 
^ 1-^0 

OJq — Ltp- 


1/2 


(19) 


The frequency dependence of the wave number together 
with those of the phase velocity and effective density are 
plotted on normalized scales in Figs. 3(d) and (e) for the 
case wi/wo = 2. In the frequency range ujq < oj < wi, 
Peff is negative and the waves are damped. Because this 
system has more degrees of freedom than the membrane 
system previously discussed, it has a more complicated 
dispersion relation. 


spring in each unit cell attached to the membrane, as 
shown in Fig. 4(a). This provides more degrees of free¬ 
dom than even the previous example of the mass-and- 
spring analogy of the solid-matrix acoustic metamate¬ 
rial. The applied force F = ApS on a unit cell from 
the pressure gradient acts on the part of the unit-cell 
mass M = poSD -(- Mmem made up as before of the sum 
of the masses of the air and the membrane. The hid¬ 
den force Fh = — v) — km^ is the combination of 

the forces from the spring —k{^ — rj), where ^ and rj are 
the displacements of the masses M and Me, respectively, 
and the force —km^ from the membrane spring constant. 
The equations of motion for sinusoidal excitation of a 
single unit cell are M^ -(- k{^ — p) + km^ = Fq exp{—iujt) 
and Me 77 + k{rj — ^) = 0. The equations of motion are 
slightly different from the previous example of the mass- 
and-spring analogy of the solid-state matrix, but reduce 
to the same form when km = 0. By elimination of the 
variable 77 one can derive Meff in the following form: 


Membrane-based metamaterial including masses and 
springs 

As a final example of the hidden-force approach, con¬ 
sider a more general case of the previously analyzed 
membrane system obtained by including a mass and a 


Meff = F/l = M (1 + , (20) 

y CJq — UJ ^ J 


where loq = y/k/Me, uJi = -\- Me/M and UJ 2 = 

y^km/M. For the special case in which the membrane 
spring constant km can be neglected, we may set uj 2 = 0. 
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UNITCELLj 



FIG. 4: (a) Schematic diagram of a ID acoustic metamaterial based on membranes, springs and masses, made up of periodically- 
spaced unit cells in a tube containing air. Normalized plots (b), (c), (d) show the frequency as a function of the density, phase 
velocity and wave number, respectively. In (c) and (d) the solid and dashed lines refer to the cases of real and imaginary values, 
respectively, for the phase velocity and wave number. Normalized parameters ui/uio = 2 and u; 2 /ti;o = 1 are chosen for these 
plots. 


This results in the simpler form 

+ (21) 

which is precisely the same as Eq. (IT^ . In this special 
case, the present system is an exact analog of the mass- 
and-spring model previously considered. The effective 
modulus is still given by Bq according Eq. ([5]), because 
the pressure-volume relation, pj = —Bo/^VjJV^ is not 
affected by the addition of the mass and spring. 

Since peff = M^ff/SD and q = uj^^peff /Bq, the dis¬ 
persion relation is given by 



where p' = M/SD. The frequency spectra for the ef¬ 
fective density and phase velocity, together with the dis¬ 
persion relation, are shown in Fig. 4(b)-(d) for the case 
loi/loq = 2 and uj 2 jijJa = 1. For this choice of parameters, 
two distinct frequency bands with negative density are 
evident. 

This mechanical model also provides the expected re¬ 
sults in the limit when either the internal-spring con¬ 
stants k/2 go to infinity or the internal mass Me is set to 


zero. In both cases the model reduces to an elementary 
mass-and-spring chain model, which, in the metamate¬ 
rial limit D <^ \, shows a constant sound velocity (i.e no 
dispersion) and positive and constant effective mass. 


To conclude this discussion of effective densities, we 
have proposed the hidden-force picture to explain why 
effective masses and densities are significantly different 
from their non-resonant average values. In this pic¬ 
ture the effective mass Mg// is obtained in terms of 
the hidden-force to applied-force ratio Fh/F as M^ff = 
M/(l -I- Fh/F). The effective mass becomes negative 
when Fh/F < —1. We demonstrated that this picture 
allows one to obtain effective masses of the unit cell and 
thereby the effective densities in a quick and easy manner 
for two established metamaterials that exhibit negative 
density as well as for a new membrane-based metamate¬ 
rial including masses and springs. 


In the next section we shall discuss with the aid of 
several examples a simple picture of how frequency- 
dependent effective moduli arise in acoustic metamate¬ 
rials. 
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EFFECTIVE MODULI: HIDDEN SOURCES OR 
HIDDEN EXPANDERS 

Here we introduce the concept of “hidden sources” of 
volume in order to understand effective modulus. This 
approach is first explained by considering the vibrational 
response of a piston connected to a chamber contain¬ 
ing either a Helmholtz resonator or a side hole. We 
then consider an acoustic metamaterial based on unidi¬ 
rectional propagation in a tube lined with Helmholtz res¬ 
onators. We go on to treat the case of a tube containing a 
combination of Helmholtz resonators and membranes—a 
generic case of a double-negative acoustic metamate¬ 
rial—followed by a similar mass-and-spring analogy. For 
mass-and-spring models we extend the concept of hidden 
sources, more appropriately termed “hidden expanders” 
of displacement in this case, by the introduction of light 
rigid trusses coupled to extra degrees of freedom for 
mechanical motion, and demonstrate an example of a 
double-negative system based on this concept. We con¬ 
clude by summarizing our approach and discussing how 
to tell at first glance what produces effective density and 
what produces effective modulus. 



FIG. 5: (a) Showing the basic concept of a hidden source, in 
which volume AVh of air is introduced into a chamber while 
depressing a piston to change the chamber volume by AU. 
Here the changes as shown correspond to negative values of 
AV and AVh- (b) shows a normalized plot of B~^j: for this 
system as a function of the hidden source of volume AVh- 

behavior is analogous to that observed in the case of hid¬ 
den forces. (The curve in Fig. 5(b) is identical in shape 
to that in Fig. 1(c).) Here one can see that a hidden 
source is represented by an introduced volume of air. 


Concept of a hidden source 


Single Helmholtz resonator 


Consider a chamber and piston containing air as well 
as a point where air can be introduced or removed. This 
point, not apparent to the operator moving the piston, 
constitutes the origin of what we call a hidden source 
or sink of volume. A schematic diagram of this setup is 
shown in Fig. 5(a), where we represent a small change 
in which the piston is displaced to perturb the chamber 
volume by AV at the same time as a volume AVh of 
air is introduced (measured at the equilibrium pressure 
before the change). The change in pressure p inside the 
chamber is given by 


^ AV + AV^ 
P = -Bo --- 


(23) 


The effective bulk modulus B^ff only depends on the 
observable volume change AV, so, in accord with the 
definition of Eq. ([ 6 |), p = —BgffAVjV defines B^ff for 
this system: 


Beff =Bo(l + . (24) 

By analogy to the definition of Mg// in Eq. ([2), B^/f 
depends on the ratio of the hidden source AV?t to a more 
easily observable quantity, here the change in chamber 
volume AV. A plot of vs AVh/AV is shown 

in Fig. 5(b). Notably, as AVi approaches — AV, B^ff be¬ 
comes zero (and compressibility B~j^ becomes infinite). 
Also, Beff becomes negative when AVh < — AV. This 


A prototype of a system that can exhibit negative Bef / 
is a piston and chamber with an attached Helmholtz res¬ 
onator of volume Vh, as shown in Fig. 6 (a). The neck of 
the resonator is assumed to have area Sr and effective 
length!^ V. This system can be regarded as a piston 
with a chamber of equilibrium volume V, with pressure 
variations p subject to a hidden source of volume that 
is governed by the pressure variation pn = BoShv/Vh 
inside the Helmholtz resonator, where rj is the displace¬ 
ment of the air plug of the resonator neck in the out¬ 
ward direction with respect to volume V. In the case 
in which the external driving force on the piston is sinu¬ 
soidal at angular frequency w, i.e. p = po exp(—iwt) and 
77 = 770 exp(— zwt), the acceleration 77 can be calculated 
from the equation of motion, poSrI'^ = Snip — Ph), 
yielding 


VhPo 1 
BoSh 1 


(25) 


where, for this case, wq = \/BqSh/{Vh pol') is the clas¬ 
sical Helmholtz resonator frequency, This treatment 
does not impose a limit on the ratio VnjV, although, 
usually, Vh < V. The volume of the neck is, however, as¬ 
sumed to be much smaller than V for the lumped-element 
treatment of the motion of the air inside it to apply. 
Higher-order resonances of the system are neglected in 
this approach. Knowing the displacement 77 allows us to 












calculate the hidden source AI 4 ,: 


AI4 = ShV = 


VhP 


Bn 


1 — 


(26) 


Using the harmonic expression AU = AVb exp(—iwt), 
Eq. (I^5)l . and the definition of i?e// in Eq. (1^^ . we 
obtain, for the compressibility = —AVb/(Upo)i 


^e 7 / = V 


1 + 


(27) 


where wi = wo\/(l + VhIV). The form of this equation 
is identical to that for Mg// for the mass-and-spring anal¬ 
ogy of the solid-matrix metamaterial and for the mass- 
and-spring membrane-based metamaterial in Eqs. m 
and m, respectively. A plot of frequency vs B^ 


eff 


IS 


shown in Fig. 6(b) on a normalized scale for the case 
uji/uiQ = 2. If the sinusoidal pressure amplitude pn is 
assumed be the imposed quantity, the volume amplitude 
Vo is determined by B~^j:. The volume will oscillate with 
large amplitude for lj approaching ujq from below because 
B~^j: becomes very big. In contrast, at a; = wi when 
B~^j=0 the system becomes infinitely rigid and the vol¬ 
ume amplitude becomes zero. The abrupt shift of the 
phase of the volume variations by tt with respect to the 
driving pressure is also predicted as the frequency passes 
through the resonance ujq ■ Since the sign of Bgf f changes 
at oj = Wo, the region of negative B^ff between wq and 
wi exhibits wave damping. Negative B^ff, a consequence 
of AV/i > —AU, thus implies in the case of sinusoidal 
excitation that the hidden source AU/i is not only in an¬ 
tiphase with but also has a magnitude larger than that of 
AU. In the limit w = 0, B^ff = Uo(l + Vh/V)~^, which 
is reduced from the expected value Bq owing to the in¬ 
crease in the total effective volume from U to U -I- Vh- 
In the limit w = 00, 7?e// = Bq because the flow of air to 
and fro from the Helmholtz resonator is effectively frozen 
owing to the inertia of the air plug in the resonator neck. 

The case of a side hole instead of a Helmholtz res¬ 
onator, as shown in Fig. 7(a), is also one of practical 
interest in acoustic metamaterial design. P 


11, For a 


side hole one may set wq = 0 because the Helmholtz res¬ 
onator stiffness (i.e. spring constant), kn = ShBq/Vh, 
vanishes. Equation (1^51) is modified to 


Vo = 


Po 


oj'^Pol' 


giving 


and 


AVh = ShV = ~ 


uj'^pnV 


B:h = h-M 1 - ^ 

U)^ 


’eff - ^0 


(28) 


(29) 


(30) 



FIG. 6: (a) shows a prototype of a system that can exhibit 
negative B^jj, consisting of a piston and chamber with an 
attached Helmholtz resonator of volume Vh- (b) shows a nor¬ 
malized plot of frequency as a function of the compressibility 
^7/ system for the case wi/wq = 2. 




COMPRESSIBILITY 


FIG. 7: (a) shows a prototype of a system that can exhibit 
negative B^ff, consisting of a piston and chamber with a side 
hole, (b) shows a normalized plot of frequency as a function 
of the compressibility B~^j> for this system. 


where W 2 = BqSh/(VpqI')- The resonance frequency W 2 
depends on U instead of Vh in this case. Equation (l30l) 
is the exact analog of Eq. ([S]) for the case of the effective 
mass of a membrane-based metamaterial. A normalized 
plot of the frequency dependence of B~jj: is shown in Fig. 
7(b). Negative B^ff is exhibited up to w = 012 . In the 
limit w = 0 , Beff = 0 , as expected since the oscillating 
air is completely free to escape from the chamber in this 
case. In the limit w = 00, Beff = Bq because the flow 
of air to and fro from the side hole is effectively frozen 
owing to the inertia of the air plug (as was the case with 
the Helmholtz resonator). 


Helmholtz-resonator-based acoustic metamaterial 


We are now in a position to consider the example 
of a ID acoustic metamaterial based on an array of 
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Helmholtz-resonators spaced at regular intervals in a air- 
filled tube 


resonators spaced at regular intervi 

i, [ffii, In, 0,0, [ 23 , s m 


4lH48| 


shown schematically in Fig. 8(a). The unit cell con¬ 
sists of a section of cylindrical tube containing a single 
Helmholtz resonator attached to the tube wall. 

We first note that the effective density is equal to the 
density of air, Pe// = P 01 as is evident from the pre¬ 
viously treated case of a tube containing membranes. 
(As there are no membranes here, the previously treated 
case, except with = 0, applies.) Here the unit cell 
volume change Al^ can be considered to depend on 
the non-equilibrium particle displacements and 
at the two unit cell boundaries, which act like pistons: 
Al^- = —^j)- (The definition of here is distinct 

from that used for the membrane-metamaterial. Here 
it refers to the acoustic displacement at the left-hand 
boundary of the unit cell, rather than that of the center 
of mass of the cell.) However, in contrast to the situation 
for the membrane-based metamaterial [Eq. ([T])], the av¬ 
erage pressure change in the unit cell now also contains 
a contribution from the hidden source: 


Pj = 


AH, -H AH 


hj 


V 


which, in the present case, can be expressed as 
Pj = - G) + H'S'ff], 


(31) 


(32) 


where r/j is the outward displacement of the Helmholtz- 
resonator air plug. From Eq. (|24ll . or, equivalently, using 
the definition 


Pj — Beff 


SD 


we obtain 


Beff = So I 1 + 


SHVj 




(33) 


(34) 


Introducing sinusoidally-varying quantities as before, 
making use of Eqs. (ESI) . (ISTl) and (IMll . and again assum¬ 


ing that D X, one again obtains Eq. ^ for At 

the resonance frequency uj = uJi, the effective modulus 
and the phase velocity become infinite. By analogy with 
zero-density metamaterials, this situation is useful for ap- 
plications in extraordinary acoustic transmission. |3Cll l3l| 
To derive the dispersion relation in the limit S <C A, 
we make use of p = —Bi.ffdu/dx and —dp/dx = pqU 
by analo gy with Eqs. ([5]) and (fTUll to derive the wave 
equation: [54 1 


POP = Beff 


d^p 


dx'^ ’ 

yielding = ui'^po/Beff. For this example, 

/P'0 / . . 

q = UJ 



(35) 


(36) 


The frequency dependences of the compressibility and 
the phase velocity, as well as the dispersion relation, are 
shown by normalized plots in Figs. 8(b)-(d). In the 
region of negative modulus the propagation is damped. 
This behavior is analogous to that noted for negative ef¬ 
fective mass. 

The equivalent result for the dispersion relation for an 
array of side holes instead of Helmholtz resonators is 



This has exactly the same form as the dispersion for the 
membrane-based metamaterial [Eq. (I12II ]. 

Just as one can ask the question why the addition of 
membranes produce no added rigidity for the passage 
of acoustic waves, one can also ask why Helmholtz res¬ 
onators add no effective mass. In fact effective mass is 
only added, according to Eqs. © and if there is 
an extra frequency-dependent force acting on a unit cell. 
Since there are no such hidden forces in the system of 
Helmholtz resonators but only hidden sources, these res¬ 
onators do not contribute to the effective mass. Rather, 
the introduction of the Helmholtz resonators leads to an 
extra contribution —BgAVk/V to the pressure [Eq. (|32])[, 
whose effect on either side of a unit cell is the addition 
of a pair of equal and opposite forces. This pair of forces 
obviously does not contribute to an imbalance in the net 
force on a unit cell, as is required for the addition of 
effective mass [Eq. (|4|)[, but instead leads to a net com¬ 
pression or expansion of the unit cell, i.e. resulting in a 
contribution to the elastic modulus. 


Acoustic metamaterial with both Helmholtz 
resonators and membranes 


We now consider a ID acoustic metamaterial based 
on Helmholtz-resonators combined with membranes in 
an air-filled tube, [3, 0, 0, I0 [ a prototype system for 
double-negative behavior, as shown in Fig. 9(a). The 
membranes and Helmholtz resonators are positioned al¬ 
ternately inside an air-filled tube. Consider the unit cell 
j sketched in the dashed line in Fig. 9(a) that ends just 
after a membrane. This unit cell is chosen so it can ap¬ 
ply to both the analysis of the membranes and Helmholtz 
resonators. The equation 


{Pj - Pj+i)S - km^j = Miij, (38) 

valid for the previously-treated case of membranes only, 
and 

Pj = ~ + VjSh], (39) 

i.e. Eq. (|32l) . derived for Helmholtz resonators only, still 
apply, where the mass M again refers to the lumped mo¬ 
tion of the unit cell. The change in the choice of unit cell 
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FIG. 8: (a) Schematic diagram of a ID Helmholtz-resonator-based acoustic metamaterial, showing periodically-spaced res¬ 
onators in a tube. Normalized plots (b), (c), (d) for this system show the frequency as a function of the density, phase velocity 
and wave number, respectively, for the case wi/wo = 2. In (c) and (d) the solid and dashed lines refer to the cases of real and 
imaginary values, respectively, for the phase velocity and wave number. In the horizontal axes we use velocity vo = i/Bo /po • 


affects the definition of the quantity in Eq. (l38l) , which 
now refers to the acoustic displacement at the left-hand 
side of the unit cell j rather than to that of the cell cen¬ 
ter of mass. However, the difference in these definitions, 
only affecting distances ^11/2, does not lead to a change 
in the final results for effective physical properties and 
the dispersion relation. The above equations separately 
determine the effective density and modulus according to 
Eqs. ([H) and (|S]), so we arrive at expressions for peff and 
B~f j: in exactly the same form as those in Eqs. m and 
(l?7)) . respectively: 

(40) 


B 


-1 

eff 




UJ 


IH 


^OH 



(41) 


where we have added the labels m for membrane and H 
for Helmholtz resonator to remove the ambiguity in the 
definitions uJom = \/k^/M, lo qh = \/B^Sh/{VhP ol') 
and bJiH = ^oH\/ (1 + VhIV). By analogy with Eq. 
(HU, the wave equation. 


Peffi = Beff 


dx^ ’ 


(42) 


leads to the dispersion relation q = p^ff /E^ff, i.e.. 



The frequency dependence of the phase velocity, as well 
as the dispersion relation, are shown by normalized plots 
Fig. 9(d), (e) for the case ujih/ujoh = 2 and LUom/i^oH = 
3, for which a region of double-negativity exists from 
1 < oj/uJoH < 2. In this region the phase velocity is oppo¬ 
site to the group velocity. In 2D and 3D such materials 
are expected to be important in focusing applications. 
The dispersion relation for the case of side holes instead 
of Helmholtz resonators can be easily found by the use of 
Eq. (|30p for B~jj: instead of Eq. (E71) . We now turn to 
the case of a mass-and-spring model exhibiting an effec¬ 
tive modulus or exhibiting both an effective density and 
modulus. 


Mass-and-spring analogy for an acoustic 
metamaterial exhibiting negative modulus or 
double-negative behavior 

Consider a ID model consisting of masses and springs 
connected to light rigid hinged trusses coupled to extra 
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UNIT CELL] 




FIG. 9: (a) Schematic diagram of a ID acoustic metamaterial based on membranes and Helmholtz resonators, showing 

alternately-spaced elements in a tube. Normalized plots (b)-(e) for this system show the frequency as a function of the 
density, phase velocity and wave number, respectively, for the case cjirr/oioir = 2 and uJom/ijJoH = 3, for which a region of 
double-negativity exists. In (d) and (e) the solid and dashed lines refer to the cases of real and imaginary values, respectively, 
for the phase velocity and wave number. 


degrees of freedom for mechanical motion, as shown in 
Fig. 10(a). We assume all mechanical displacements are 
much smaller than the truss lengths. This type of hinged 
truss system was previously proposed together with ex¬ 
tra springs to generate a negative modulus. (4fl| but we 
have simplified the model to a convenient bare minimum 
here. Springs of constant 2fco connect a mass M in the 
unit cell to the truss systems. The square truss system, 
consisting of four members, is connected above and below 
to straight trusses, which in turn are connected to two 
masses m that are constrained (by rails) to move only 
in the vertical direction. These ideal (i.e. massless and 


frictionless) hinged trusses ensure 1) that the same mag¬ 
nitude of force is exerted on the springs either side, and 
2 ) that the vertical displacements y of the masses m are 
exactly mirrored by the horizontal displacements of the 
sides of the trusses attached to the springs. The presence 
of the two vertically-oriented straight trusses allows the 
square truss system to be free to move horizontally. The 
springs provide the applied force F = Fj — F)+i on a 
particular unit cell j, where, from Eq. a, 

Meff = (44) 

si 
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However, from the force transmission properties of the 
truss system, that ensure that the compressive forces in 
the springs on either side of it are equal [see Fig. 10(a)], it 
is clear that the acceleration of mass M is simply given by 
^3 — i ^3 ~ ^ 3 +i)/^- Therefore, from Eq. (Hill . Meff = 
M for this model. 

To find the effective modulus, first consider the force 
Fj on the left-hand side of unit cell j. This depends 
on the displacement fij-i — yj-i of the right-hand side 
of the square truss of unit cell j — 1 as well as on the 
displacement of mass M in unit cell j, where yj and 
—yj are the displacements of the top and bottom masses 
m in unit cell j, respectively, and y,j is the horizontal 
displacement of the center of mass of the square portion 
of the truss system: 

Fj = 2fco(Mj-i - 2/1-1 - (45) 

The compressional force F)_|_i in the spring to the right 
of mass M in unit cell j [see Fig. 10(a)[ is given by 


that shown in Fig. 7(b). So the model of Fig. 10(a) is the 
mechanical analog of an air-hlled tube with periodically 
arranged side holes. At w = 0 the effective modulus 
is zero because the truss system provides zero effective 
spring constant in this limit. In contrast to the mass- 
and-spring model of Fig. 3(b), the springs in the model 
of Fig. 10(a) cannot support any tension or compression 
in their equilibrium position, i.e. the springs should work 
in both tension and compression in the present case.fh^ 
At w = oo the effective modulus becomes equal to ko/D, 
identical to that of the previously considered mass-and- 
spring model, because in this limit the masses m do not 
move and the square truss plays the role of a massless, 
rigid connector, resulting in two springs of constant 2kQ 
in series that are equivalent to a single spring of constant 
ko- In this limit the effective modulus is simply ko/D. 

The question arises of how to interpret this mechanical 
model in terms of the hidden-source picture. The analo¬ 
gous equation to Eq. (l23|l for the mass-and-spring model 
is, by comparison with the definition of Eq. (H, 


Fj +1 = 2fco(0 - - y 3 )■ (46) 


Decreasing the indices by 1 yields an alternative expres¬ 
sion for Fj: 


Fj = 2fco(^j-i - pj-i - yj-i). (47) 


Comparing Eqs. (PSI) and ([4711 . we obtain 


Pj = 


^3 + 0-1 

2 


(48) 


The massless square truss system moves by this amount 
to ensure equal displacements in the springs on either 
side of it and thus maintain the force balance on it. (Any 
massless system must by definition have a net zero force 
on it to avoid an infinite acceleration.) Eliminating fij 
from Eq. ([45|l 


Fj — ko{^j-i ^j 2yj-i). (49) 

From the properties of the truss system, force Fj is trans¬ 
mitted to the masses m. For sinusoidal motion at angular 
frequency w, Fj = —moj^yj-i. This allows yj-i to be ex¬ 
pressed as 


2/i-i = - 


^3 0-1 

2[\-oFlujir 


(50) 


where U3i = ^/Fk^jm, for this case. The definition of the 
effective Young’s modulus, Eq. (HU, then leads to 


F=-D‘^ — ^—— ^ (52) 
^ D D D D ' ^ ^ 

where k^/D is the modulus in the absence of the truss 
system and (^j — ^j-i)h is an extra displacement that 
we term a “hidden expansion”. By comparing Eq. (l52l) 
with Eqs. ([45l) and (l50|l . one finds that {^j — ^j-i)h = 
2yj-i. In other words, the hidden expansion is precisely 
equal to the extra horizontal displacement introduced by 
the truss system. The mass-and-spring analogy of the 
hidden source concept is thus, quite naturally, a hidden 
expander. 

The acoustic dispersion relation can be derived by 
analogy with the previously considered mass-and-spring 
model. 


= 2ko{y,j-i — yj-i — ^j) — 2ko{^j — y,j — yj), (53) 
Provided that D <C A, this equation reduces to 

= (54) 

Compared to the previous case of Eq. (HU), there is an 
extra term in dy/dx owing to the truss system. A differ¬ 
ential equation involving y can be derived by noting that 
m{yj — jjj-i) = Fj+i — Fj, making use of Eqs. (1451) . ([46l) 
and (I48L assuming D A, and then integrating once 
over the coordinate x: 

my = -koD- - 2koy. (55) 


’-1 

ko 



(51) 


The effective modulus varies with frequency in exactly 
the same way as a tube containing an array of side holes 
[as in Eq. (HOI) ]. The frequency variation is the same as 


Substituting parameters with temporal variations ac¬ 
cording to exp [/(gs—tut)] as before leads to the dispersion 
relation q = ujy/p/Eeff, where p = M/D^: 


q = Lo 




1/2 


(56) 
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UNIT CELLj 


FIG. 10: (a) Mass-and-spring model that can exhibit negative modulus, realised by periodically-arranged masses and springs 
connected to light rigid trusses coupled to extra degrees of freedom for mechanical motion, (b) Mass-and-spring model that 
can exhibit double-negative behavior. Fj and Fj+i are the compressional forces in the springs on the left-hand side of and to 
the right of unit cell j, respectively. The ideal light rigid truss system containing two attached masses m is hinged in such a 
way as to ensure that the compressional (or tensional) forces in the springs on either side of it are equal. The masses m are 
constrained to move along fixed vertical rails. 


which has exactly the same form as Eqs. (TT^ and ((571) . 

Let us now turn to a more general ID mass-and-spring 
model that can exhibit double-negative behavior, as il¬ 
lustrated in Fig. 10(b). We have combined the model of 
Fig. 10(a) with that of Fig. 3(b). The analysis proceeds 
in exactly the same way as for these two cases: Meff and 
Eeff are given by Fqs. (fT3)) and (IhTl) . respectively, and 
the dispersion relation, q = uj yjpef / /Eef /, becomes 


q = 



2 2 
^ 1-^0 

Wq — 




Somewhat coincidentally, this has precisely the same 
form as that of Fq. (H5)) for the case of an air-filled tube 
containing a periodic array of membranes and Helmholtz 
resonators. The frequency dependence of the term aris¬ 
ing from the hidden expanders is the same as that for the 
membranes (which give rise to hidden forces), whereas 
the frequency dependence of the term arising from the 
hidden forces has the same form as that for the Helmholtz 
resonators (which give rise to hidden sources). The fre¬ 
quency spectra of Pe// and Fg// are analogous to the 
plots of Fig. 9(b), (c), the only difference being that the 
roles of the effective mass and modulus are reversed. The 
plots for jupl and |g| are the same as for Fig. 9(d) and 
(e) for equivalent dimensionless parameters. 


We conclude this section by emphasizing that effective 
mass and modulus are not just theoretical constructs but 
also experimentally measurable quantities. According to 
their definitions in Fqs. o, (ED and (USD, it suffices in 
principle to put pressure, force or displacement sensors 
at the appropriate points inside the metamaterial, and 
then the effective parameters can be experimentally de¬ 
rived. In contrast, merely measuring the dispersion re¬ 
lation, that depends on the combination p^ff/Beff or 
Peff/Egff, will not in general be sufficient to distinguish 
effective mass from effective modulus. 


CONCLUSIONS 

In conclusion, we have proposed the concepts of hid¬ 
den force and hidden source of volume to respectively 
account for the effective densities and moduli of acous¬ 
tic metamaterials. The superficially strange concepts of 
negative density and modulus are naturally accounted 
for in this picture when the hidden force/source oper¬ 
ates in antiphase to and is bigger in magnitude than the 
force/volume-change engendered in its absence. We illus¬ 
trate our approach in ID for well-known air-based meta¬ 
materials involving membranes, Helmholtz resonators or 
side holes with the inclusion of the new case of an array of 
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masses attached to membranes by springs. We also intro¬ 
duce examples based on generic mass-and-spring models, 
in which case the concept of a hidden source of volume 
is replaced by the concept of a hidden expander of dis¬ 
placement. 


Deciding at first glance what contributes to an effective 
density or to an effective modulus depends on the sys¬ 
tem. For air-based acoustic metamaterials, membranes 
only contribute to the effective density whether or not 
they have an attached mass and spring. The reason 
for this is that they only involve local, time-dependent 
hidden forces that provide a net force on the unit cell. 
A similar result applies to the mass-and-spring model 
representing a solid-matrix acoustic metamaterial based 
on heavy spheres in a soft matrix. In air-based acous¬ 
tic metamaterials the effective modulus arises because of 
time-dependent hidden sources of air volume associated 
with Helmholtz resonators or side holes, that only pro¬ 
duce pairs of forces acting equally and oppositely on ei¬ 
ther side of the unit cell (i.e. resulting in a zero net force 
on the unit cell). A similar result applies to the mass- 
and-spring models containing light rigid hinged trusses 
attached to masses. Although not discussed here, the in¬ 
clusion of extra degrees of freedom in mass-and-spring or 
solid-state models, e.g. through rotations, can also result 
in an effective modulus through the production of pairs 
of equal and opposite forces. [ISl l49|, 


The extension of these ideas to 2D or 3D should be 
straightforward, whether for fluid, solid or multi-phase 
metamaterials. [13| In these cases, due regard should be 
taken of possible anisotropic properties and of the dif¬ 
ferent acoustic modes of propagation such as longitu¬ 
dinal, shear or flexural, for example. For simplicity in 
our treatment we have restricted our attention to fun¬ 
damental resonances, namely of membranes, Helmholtz 
resonators or side holes, although the inclusion of higher- 
order resonances in the general framework presented is 
possible. We have also ignored material damping. This 
has the advantage of highlighting the damping caused 
by the intrinsic metamaterial properties. Frequency re¬ 
gions for single-negative-parameter behavior exhibiting 
damping can be referred to as band gaps, although they 
should not to be confused with those arising from purely 
phononic effects (generally observed at higher frequen¬ 
cies). Frequency regions for double-negative behavior are 
of particular interest because of their potential for high- 
resolution focusing. Finally, acoustic metamaterials will 
probably be entering the application stage as commer¬ 
cial products in sound control in the near future, and we 
hope that this paper will aid in accelerating progress in 
this regard. 
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